Damage and recovery of fiber Bragg grating under radiation environment*

Project supported by the Project for the State Key Laboratory of Optoelectronic Materials and Technologies of China (Grant No. 09010-32031708) and the Project for Zhuhai Key Laboratory of Center for Space Technology of China (Grant No. 71000-42080001).

Wen Shi-Zhe1, 6, Xiong Wei-Chen3, 6, Huang Li-Ping4, 6, Wang Zhen-Rui4, 6, Zhang Xing-Bin5, 6, He Zhen-Hui2, 5, 6, †
School of Material Science and Engineering, Sun Yat-Sen University, Guangzhou 510006, China
State Key Laboratory of Optoelectronic Materials and Technologies, Sun Yat-Sen University, Guangzhou 510275, China
School of Engineering, Sun Yat-Sen University, Guangzhou 510006, China
School of Physics, Sun Yat-Sen University, Guangzhou 510275, China
School of Physics and Astronomy, Sun Yat-Sen University, Zhuhai 519000, China
Center for Space Technology, Sun-Yat Sen University, Zhuhai 519000, China

 

† Corresponding author. E-mail: stshzh@mail.sysu.edu.cn

Project supported by the Project for the State Key Laboratory of Optoelectronic Materials and Technologies of China (Grant No. 09010-32031708) and the Project for Zhuhai Key Laboratory of Center for Space Technology of China (Grant No. 71000-42080001).

Abstract

To develop the application of fiber Bragg gratings as temperature and strain sensors in space environments, it is necessary to understand the effect of high-energy radiation on the performance of the fiber Bragg grating. We performed an experiment involving Co60-γ ionizing irradiation with a total dose of 1.01 × 106 rad on two Ge-doped single-mode fiber Bragg gratings with central wavelengths of 825 and 835 nm, respectively. We found that, with the increase of radiation dose, the redshift of the peak wavelength of the reflection spectrum of the fiber Bragg gratings indicated the increase of the refractive index and the number of color centers. After irradiation, the refractive index decreased with the decreasing number of color centers. We analyzed the influence of ionizing irradiation on the transmission performance of the fiber Bragg gratings using a color-center model, which explained the experimental results. The proposed model was used to determine the creation rate and annihilation rates of the color center, which are foundational data for using the fiber Bragg gratings in space applications.

1. Introduction

Fiber Bragg gratings are regarded as filters or mirrors for lasers with the corresponding wavelength. Through internal writing and ultraviolet phase mask exposure, the refractive index of a grating exhibits a periodic distribution. With the above features, Bragg gratings can be used as optical devices, having the advantages of small size, lightness, high strength, anti-electromagnetic interference, corrosion resistance, easy coupling with fibers, small coupling loss, and low cost. In satellite laser communications, fiber Bragg gratings can be used as a narrow-band filter, for wavelength-division multiplexing (WDM) and as a demultiplexing unit of WDM systems, and as an optical branching multiplexer.[13] Additionally, fiber-optic sensors can detect the working status of satellites or nuclear facilities through the displacement, speed, temperature, and other measurements of the reaction of the facilities to ensure their operating health.[46] Fiber Bragg gratings have better multiplexing performance and can be used for sensing and transferring data, whereby it is easy to form a sensor network.[7]

To evaluate the performance of fiber Bragg gratings in high-energy radiation environments, e.g., in space or nuclear facilities, experiments should be performed in similar environments to understand the effect of the radiation on the grating performance. Because of high cost, long period, and difficulty, a Co60-γ source is typically used to simulate the radiation environment. Ferdinand[4] and Gusarov[5,6] conducted γ and neutron radiation experiments on a fiber Bragg grating on a germanium fiber and observed the increase of the peak wavelength and the refractive index of the light reflected from the grating. The phenomenon was attributed to the color centers created by ionizing radiation. Finally, they did not conduct an in-depth analysis of the regularity and tendency of the color centers.[46]

A common feature of aerospace and nuclear reactors is the harsh radiation environment, which can damage the performance of fiber-optic components.[1,46,8,9] Radiation is deposited on the material mainly in two ways: the displacement of atoms and ionization. Owing to the strong covalent-bond structure in silica optical fiber, it is very difficult for radiation particles to displace the atoms. Therefore, the ionization effect is the major cause of defects.[10] For silica fiber, the ionization effect shifts an electron to the neighboring atoms and thus creates a defect called a color centers. The created color centers change the optical absorption spectrum,[1115] modifying the refractive index of the fiber Bragg grating.[1,6,16]

In recent years, among several studies investigating the space radiation effects on optical fibers,[1719] few paid attention to the radiation damage to fiber Bragg gratings, and none investigated the recovery from such damage.

In this study, to examine the radiation effects on the fiber Bragg grating, an experiment involving γ60Co irradiation, along with post-irradiation annealing at room temperature, was performed. According to the model for the creation and annihilation of color centers, the creation rate and the annihilation rate of color centers were fitted from the experiment data, and the shift of the central wavelength caused by the irradiation was predicted and set as the systematic error.

2. Radiation-created color-center model

Ma et al.[16] introduced a model of color-center creation and degradation, where the created color centers can be degraded back to the intrinsic defects by the radiation, and the degradation rate of the color center is proportional to the radiation rate. However, this assumption was not supported by their experimental results, which indicated that the properties of the fiber recovered gradually after radiation, indicating that the color-center concentration decreased when there was no radiation.

We consider a different model of color-center creation and annihilation. Color centers can be created by high-energy radiation but are unstable and can be annihilated by thermal annealing. That is, under constant irradiation over a long period of time, the creation rate becomes equal to the annihilation rate, and the color-center concentration becomes constant. Consider the combined effect of the color-center creation rate kc, which is proportional to the radiation dose rate Rd and to the concentration of atoms in regions without color center N0, as well as the color-center annihilation rate ka, which is proportional to the available concentration of color center n(t). The variation of the color-center concentration in the fiber can be expressed as

where κc is a coefficient of the dose rate, and t is the irradiation time. kc and ka can be obtained by fitting the measured data.

When kc ≠ 0, the solution of equation (1) gives

By replacing Rdt in Eq. (3) with the irradiation dose D = Rdt, equation (3) can be rewritten as

When irradiation is stopped, kc = 0, and equation (1) becomes

with the solution
where n(t0) is the color-center concentration immediately after the irradiation is stopped.

In Smakula’s equation,[20] the absorption coefficient is proportional to the color-center concentration

According to the Kramers–Kroning causality law,[21] the relation between the absorption coefficient and the refractive index can be expressed as
where λi is the center wavelength of the i-th absorption band, Δαmax is the variation of the peak absorption, Δλ1/2 is the full width at half maximum (FWHM) of the absorption band, and λ′ is the wavelength of interest.

Consider the simplest case: for an optical fiber Bragg grating around 800 nm, the ionizing effect of electron radiation on the grating mainly produces one type of color center,[22] which leads to an absorption spectrum of Gaussian approximation. When the FWHM of the absorption band is a certain value, from Eqs. (4), (8), and (9), we obtain the following relationship for the fiber Bragg grating under irradiation:

where η is a constant associated with the type of color center and the absorption spectrum, and φ is a constant associated with the initial color-center concentration and η.

From the Bragg grating equation

we have
where Δλmax is the peak wavelength variation, Λ is the grating period length, and λB is the Bragg wavelength.

After radiation, Δnz can be obtained from the measured Δλmax, and then η can be obtained by fitting Eq. (10). Finally, the effect of high-energy radiation on the fiber Bragg grating sensor can be predicted.

3. Experiments
3.1. Apparatus and method

The fiber Bragg gratings used in this experiment were procured from FBGS. The gratings were inscribed on a germanium-doped (10%–20%) single-mode fiber via the phase mask method. The fiber core diameter and cladding diameter were 4.5 and 125 μm, respectively, and the coating material was modified organic ceramics, with the temperature sensitivity of the Bragg wavelength of the grating being 5.5 pm/°C (according to the data sheet for the fiber Bragg grating). The parameters of the grating sample are shown in Table 1, and the grating structure under the microscope is shown in Fig. 1.

Fig. 1. (color online) Fiber Bragg grating structure.
Table 1.

Parameters of optical fiber Bragg gratings.

.

The radiation source used in the experiment was the 60Co-γ radiation source from Guangzhou Huada Biological Co., Ltd. The grating section of the optical fiber was exposed in a 60Co-γ-ray radiation chamber, while the other part of the optical fiber was protected by a 3-mm-thick lead tube. Near the grating section, a low-range potassium dichromate dosimeter was used to measure the radiation dose. The irradiation period of 2 h for each data point was not strictly controlled, with an uncertainty of 10 min. After each irradiation, the reflection spectrum of the fiber sample was measured. The interval between two irradiations was hours on the same day, or overnight, depending on the working time of the company. The total radiation dose was 1.01 × 106 rad. A K-type thermocouple was used to monitor the ambient temperature.

A sketch of the experimental setup is shown in Fig. 2. The built-in wavelength of the light source ranges from 800 to 860 nm, and the intensity is approximately 0.3 mW. The emitted light passes through the fiber optic splitter and then is split into two beams of equal intensity to the two fiber Bragg grating. The two specific wavelengths of light are reflected into the spectrometer. Through an array of charge-coupled devices (CCDs), the reflected light intensity at its wavelength is transformed into an electrical signal to record the reflection spectrum.

Fig. 2. (color online) Experimental setup for the radiation experiments of fiber Bragg gratings.
3.2. Radiation dose error analysis

The main sources of measurement uncertainty are the magnitude of the change of the light absorption before and after irradiation, ΔA, and the ambient temperature T. The sum of squared uncertainties is multiplied by the dose response conversion factor K (= 8.33) to obtain the total uncertainty of the measurement. The K-type thermocouple temperature measurement accuracy is ±0.5 K, and the absorbance accuracy A is 0.001; the dose uncertainty can be calculated using the following equation:

The relative error of the calculated dose is

Table 2 shows the absorbed dose errors at different dose rates. The relative error varies from 0.39% to 0.77%.

Table 2.

Absorbed dose error.

.
4. Results and discussion
4.1. Effect of radiation dose on reflection spectrum

Tables 3 and 4 show the variation of the peak refractive wavelengths of the 825- and 835-nm fiber Bragg gratings, respectively, according to the radiation dose. As the temperature of the fiber Bragg grating changes by 1.0 °C, the peak wavelength changes by 5.5 pm. Considering the effect of the ambient temperature, the peak wavelengths of the fiber Bragg gratings were corrected. Figures 3 and 4 show the reflection spectrum in two grating samples before and after irradiation, where a redshift of the peak wavelength is obviously observed. The 825- and 835-nm peak wavelengths increase by 0.0706 and 0.0742 nm, respectively. The FWHM of the reflected wave does not change significantly with the increase of the irradiation dose (see Tables 5 and 6). With the increase of the irradiation dose, the intensity of the reflected wave of Bragg grating first increases and then decreases (see Tables 5 and 6), indicating that the light-transmission characteristics of the fiber were impaired by the radiation.

Fig. 3. (color online) Reflective spectrum variation of the 825-nm optical Bragg grating pre-radiation and post-radiation.
Fig. 4. (color online) Reflective spectrum variation of the 835-nm optical fiber pre-radiation and post-radiation.
Table 3.

825-nm Bragg grating wavelength variation with respect to radiation dose.

.
Table 4.

835-nm Bragg grating wavelength variation with respect to radiation dose.

.
Table 5.

825-nm Bragg grating FWHM and relative intensity variation with respect to radiation dose.

.
Table 6.

835-nm Bragg grating FWHM and relative intensity variation with respect to radiation dose.

.
4.2. Effect of radiation dose on refractive index of fiber Bragg grating

According to the experimental results, the difference between the peak wavelength of fiber Bragg grating before and after radiation could be obtained at different radiation doses.

The measurement accuracy of the peak wavelength of the Bragg grating reflection was 0.1 pm. Because the ambient temperature influenced the peak wavelength of the reflection, in the data-processing operation, the peak wavelength was converted into the corresponding wavelength at 20 °C. The thermocouple temperature measurement accuracy was 0.5 K. The conversion equation is

or
where λmax is the peak wavelength corresponding to the ambient temperature of 20 °C, λ is the peak wavelength actually measured, and T is the ambient temperature. Therefore, the uncertainty of λmax is

From Eq. (9), the uncertainty of the refractive-index variation is

The variation of the refractive index of the fiber Bragg grating at different doses is shown in Tables 7 and 8. The refractive index varies exponentially with time, in agreement with Eq. (10).

Table 7.

825-nm Bragg grating refractive-index variation with respect to radiation dose.

.
Table 8.

835-nm Bragg grating refractive-index variation with respect to radiation dose.

.
4.3. Color-center annihilation

After radiation (Rd = 0 at this time), the peak wavelength of the grating reflection had a recovery or annealing effect (color-center annihilation); thus, ka can be obtained by fitting the recovery-effect curve. Then, we can obtain the particular parameters in Eq. (7) with the annealing results.

After ionizing radiation at 1.0 × 106 rad, the two fiber Bragg gratings were placed at room temperature, and their refractive-index variation with time was measured. The experimental results are shown in Tables 9 and 10 and Fig. 5. The refractive index of both fibers gradually decreases with time.

Fig. 5. (color online) Refractive-index recovery with time, after radiation. The black squares and triangles represent the measured refractive-index data for 825 and 835 nm, respectively. The black line and dashed line are the curves fitted with the color-center annihilation model for 825 and 835 nm, respectively.
Table 9.

825-nm Bragg grating recovery effect.

.
Table 10.

835-nm Bragg grating recovery effect.

.

The color center created by the radiation is unstable and will gradually be annihilated. According to Eq. (11) and the experimental data presented in Tables 9 and 10, the fitted equations for the fiber Bragg gratings with Bragg wavelengths of 825 and 835 nm are y = 8.0 × 10−5 e−0.054t and y = 8.7 × 10−5 e−0.049t, respectively, where t represents time (see Fig. 5). Therefore, we obtain ka_825 = 0.054 / 24 / 60 min = 3.8 × 10−5 min−1 and ka_835 = 0.049 / 24 / 60 min = 3.4 × 10−5 min−1, respectively.

To evaluate the reliability of the recovery curves, the error of the fitted curve and the experimental data points were determined.

where Xi is the deviation of the measured data from the fitted value yi, and σ is the standard deviation of Xi. σ is 4.18 × 10−6 and 4.22 × 10−6 for 825-nm and 835-nm gratings, respectively. Therefore, ka_825 and ka_835 with random error are expressed as follows:

From Fig. 6, the deviations of the fitted data are within 3σ. According to the 3σ criterion, the fitted curve is reliable.

Fig. 6. (color online) 3σ judgment for the difference of the experimental data points from the fitted values. The black squares and triangles represent the difference between the measured refractive-index variances and the fitted recovery curve for 825-nm and 835-nm gratings, respectively. They are within the region of 3σ judgment between the black dotted lines and dashed lines for 825 and 835 nm, respectively.
4.4. Refractive-index variation with correction

According to the color-center annihilation mechanism, the refractive index recovers in each irradiation interval. Therefore, the refractive-index variations were corrected according to the fitted annihilation rate. Owing to the time interval between the placement in the irradiation chamber and the start time of the measurement, the correction time of the grating has an uncertainty of 10 min, which yields a refractive-index variation of 3.0 × 10−8. Tables 11 and 12 show the time interval and corrected refractive-index variation time. The fitted curves are shown in Fig. 7.

Fig. 7. (color online) Refractive-index variation with respect to the radiation time after recovery correction. The black squares and triangles represent the refractive-index data with respect to time (different doses) for 825- and 825-nm gratings, respectively. The black line and dashed line are the curves fitted using the color-center creation model for 825 and 835 nm, respectively.
Table 11.

Corrected 825-nm Bragg grating refractive-index variation with respect to radiation dose.

.
Table 12.

Corrected 835-nm Bragg grating refractive-index variation with radiation dose.

.

According to ka from the annealing experiment and the fitted curve in Fig. 7, the relationship between the variations of the effective refractive index of the optical fiber with 825- and 835-nm Bragg gratings and their absorbed dose are as follows:

From Eqs. (10) and (22)–(25), ka_825 + κc_825Rd = 1.9 × 10−3, and ka_825 + κc_835Rd = 2.4 × 10−3. As Rd = 1400 rad/min, ka_825 = 3.8 × 10−5 min−1, ka_825 = 3.4 × 10−5 min−1, κc_825 = 1.3 × 10−6, and κc_835 = 1.7 × 10−6.

The relationships between the refractive-index variation for Bragg wavelengths of 825 and 835 nm and the radiation time and dose rate are as follows:

4.5. Prediction of space radiation effects on fiber Bragg grating

According to the refractive-index variation and damage recovery, the 825- and 835-nm fiber grating color-center creation rates were kc_825 = κc_825Rd = 1.3 × 10−6Rd and kc_835 = κc_835Rd = 1.7 × 10−6Rd, respectively, and the annihilation rates were ka_825 = 3.8 × 10−5 min−1 and ka_825 = 3.4 × 10−5 min−1, respectively. Therefore, the average color-center creation rate was

and the average color-center annihilation rate was
Combining Eqs. (26) and (27) by averaging the related parameters gives the variation of the refractive index of the fiber Bragg grating with respect to the radiation dose and time

With Eq. (28), the application performance of the fiber Bragg gratings in radiation environments can be predicted. For example, in the case of temperature measurement, the temperature drift caused by the radiation can be predicted in low Earth orbit and Earth synchronous orbit. For a wavelength variation of 5.5 pm/°C, according to Eq. (13), the temperature drift caused by radiation is

In low Earth orbit, the dose rate is rad/min, and the temperature drift due to radiation is

According to Eq. (30), in low Earth orbit, after 125 d (180,000 min), the maximum drift of the temperature measurement is 0.018 K, which can be regarded as systematic error. Such uncertainty caused by radiation can satisfy most of the requirements for the measurement accuracy in applications.

In Earth synchronous orbit, the dose rate is Rd2 = 5.2 × 105 rad/year = 9.9 × 10−1 rad/min. Then, the temperature drift due to radiation is

According to Eq. (31), in low Earth orbit, after 125 d (180,000 min), the maximum drift of the temperature measurement is 0.71 K.

5. Summary

We built a model of color-center creation by ionizing radiation with annihilation to study the shift and recovery of the refractive index of fiber Bragg gratings made of silica under and after 60Co γ radiation. By fitting the experimental data of the refractive-index shift, we obtained the model parameters. The effect of the radiation on the performance of the silica fiber Bragg gratings can thus be quantitively predicted for applications, e.g., temperature measurement in low Earth orbit and Earth synchronous orbit. The deviation of the refractive index of the grating in a known radiation environment can be regarded as a systematic error. This makes it possible to apply the silica fiber Bragg gratings as sensors in high-energy radiation environments, e.g., in space and nuclear facilities. However, the model provides a mechanism for the creation and annihilation of color centers, which elucidates the effect of radiation on the silica fiber Bragg gratings.

Reference
[1] Gusarov A I Doyle D B Karafolas N Berghmans F 2000 Int. Symp. Opt. Sci. Technol. 253
[2] Marshall P W Dale C J Burke E A 1992 Nucl. Sci. IEEE Trans. 39 1982
[3] Ott M Plante J Shaw J Shaw J Garrisondarrin M Ott M Plante J Garrisondarrin M 1997 Aiaa Journal
[4] Ferdin P Magne S Marty V Georges T 1994 Proc. SPIE 2425 11
[5] Gusarov A Fernandez A F Vasiliev S Medvedkov O Blondel M Berghmans F 2002 Nucl. Instrum. & Methods Phys. Res. 187 79
[6] Gusarov A I Berghmans F Deparis O Fernandez A F 1999 Photon. Technol. Lett. IEEE 11 1159
[7] Marshall P W Dale C J Carts M A Label K A 1994 Nucl. Sci. IEEE Trans. 41 1958
[8] Girard S Tortech B Regnier E Uffelen M V Gusarov A Ouerdane Y Baggio J Paillet P Ferlet-Cavrois V Boukenter A 2007 IEEE Trans. Nucl. Sci. 54 2426
[9] Zhou C M Zhang F Ding L Jiang D S 2011 Laser & Optoelectronics Prog. 48 9
[10] Holmes Siedle A Adams L 2001 Oxford University Press
[11] Jiang H Chen B X Fu C S Sui G R Ji S 2010 Acta Phys. Sin. 59 7782 in Chinese
[12] Neustruev V B 1994 J. Phys.: Condens. Matter 6 6901
[13] Regnier E Flammer I Girard S Gooijer F Achten F Kuyt G 2007 IEEE Trans. Nucl. Sci. 54 1115
[14] Russell P S PoyntzWright L J Hand D P 1991 San Jose-DL tentative
[15] Xiao Z Y Luo W Y Wang T Y 2007 Acta Phys. Sin. 56 2731 in Chinese
[16] Ma J Che C Yu S Y Tan L Y Zhou Y P Wang J 2012 Acta Phys. Sin. 204 in Chinese
[17] Fox B P Simmons-Potter K Kliner D A V Moore S W 2013 J. Non-Cryst. Solids 378 79
[18] Mady F Benabdesselam M Mebrouk Y Dussardier B 2010 Ytterbium-doped Silica Optical Fibers
[19] Jin J Li Y Zhang Z C Wu C X Song N F 2016 Chin. Phys. 25 84213
[20] Gan F ZhN G Q 1989 Crystal Color Center Physics Shanghai Jiao Tong University Press
[21] H D P Russell P S J 1990 Opt. Lett. 15 102
[22] Song J M Guo J H Wang Q X Hu S L 2012 Laser & Optoelectronics Prog. 58–64